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NONSTATIONARY FLOW ABOUT A WING-AILERON-TAB COMBINATION 

INCLUDING AERODYNAMIC BALANCE 

By Theodore Tkbodohsbn and I. E. Gahsick 



SUMMARY 

The present paper presents a continuation of the work 
published in Report No. 496. The results of that paper 
have been extended to include the effect of aerodynamic 
balance and the effect of a tab added to the aileron. The 
aerodynamic coefficients are presented in a form conven- 
ient for application to the flutter problem. 

INTRODUCTION 

It is the object of this paper to present theoretical 
expressions for the forces and the moments in a uniform 
horizontal air stream on a plane airfoil performing 
small sinusoidal motions in several degrees of freedom: 
vertical motion, torsional movement about an arbitrary 
spanwise axis, aileron movement about a hinge axis 
not necessarily located at the leading edge of the aileron, 
and tab movement similar to the aileron movement. 
The solution of this problem has direct application to 
the larger problem of flutter involving these various 
degrees of freedom and, in particular, to flutter of 
tails with control surfaces, including servocontrols. 

The development of the theory is analogous with 
that of Theodorsen (reference 1) who treats explicitly 
the case of three degrees of freedom: vertical motion, 
torsional movement about an arbitrary spanwise axis, 
and an aileron movement about a_hinge axis located 
at the leading edge of the aileron. 

Since this work was originally begun, there have ap- 
peared two German papers, one by Kussner and Schwarz 
(reference 2) and one by Dietze (reference 3), that bear 
directly on the problem. A comparison of the results 
of this paper with the results of Kussner and Schwarz, 
obtained by a different development, is given in 
appendix A. 

AIR FORCES AND MOMENTS 
Figure 1 represents a wing section with two hinges, 
an aileron (rudder) hinge at x=e and a tab hinge at 



x=f. The leading edge of the wing is at x= — 1 and 
the trailing edge at x=l. The leading edge of the 
aileron is at x=c and the distance from the hinge to the 
aileron leading edge e—c is denoted by I. The leading 
edge of the tab is at x=d and the distance from the 
tab hinge to the tab leading edge/— d is denoted by m. 
The wing is undergoing the following motions with 
small amplitudes: a displacement h (velocity K) in a 
vertical direction downward; a turning about x=a, the 
instantaneous angle of attack being a; a rotation of 
the aileron about x=e, the angle of the aileron or 
rudder being jS measured with respect to the wing; 
and a rotation of the tab about x=f, the angle of the 
tab being y measured with respect to the aileron. The 



L.E. 




Fiouai l .— Representation of wlnj section with aileron and tab snowing main 
parameters. 

actual chord is considered to be of length 26, so that b 
is used throughout the analysis as a reference length. 

The procedure and method follow those of reference 1 . 
In order to avoid needless repetition of certain expres- 
sions contained in this reference, the following notation 
is frequently used. A symbol or equation followed by 
(reference 1) denotes the corresponding expression of 
reference 1. The final results will, however, be 
explicitly given independently of reference 1. 

The air forces and the moments are treated in two 
groups: the noncirculatory and the circulatory. The ex- 
pressions for the noncirculatory part consist of apparent- 
mass terms, which do not depend on the vorticity in 
the wake. The circulatory part takes account of the 
vorticity in the wake generated at the trailing edge. 
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Noncirculatory forces and moments. — The noncir- 
culatory velocity potential at the surface, associated 
with the various motions of the airfoil, is 

*=«i+*.+*i + ^+^ + *r + *i (i) 
where only 4>e, 4>p, i>t, and 4>i need be given here. 1 



<t>fi=<j)fi (reference l)+v$b- log N(x,c) 
<t>k—<t>P (reference 1) 

-i8i^[VT II ?cos- I c-(a;-c) log iVfoc)] 



(2) 



(3) 



and 



(4) 



4>'i=4>p 

with /S, e, and i replaced by and m, respectively. The 
two extra terms appearing in fa and 4>b> which contain 
the coefficient I, arise from the aerodynamic-balance 
effect, that is, the offset of the aileron hinge position 
from the aileron leading edge. The derivation is as 
follows: 

The motion of the aileron around x=e is considered 
separated into two parts, a turning around the aileron 
leading edge plus a vertical displacement of the aileron 
relative to the wing (fig. 2). The amplitude of the 





-/ 
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Figubi 2.— Representation of the motion of the aileron around x-< as separated 
Into a turning around i=e plus a vertical displacement. 

first type of motion is {3 and of the second type of 
motion is bH=b(i(e—c)=b{il. The additional potentials 
referred to are then due to the effect of the vertical 
displacement MJand the associated vertical velocity bit. 

The potential associated with bH is determined from 
a limiting case (as c'-»c) of a shape (fig. 3) located on 
the x axis from -1 to x=c and at ordinate II from 
x=c' to 3=1 and joined by a straight-line segment 
from x—c to x=c'. This potential is associated with a 
vertical-velocity distribution: 



> The expressions for the * values are to be understood prefixed by ± signs: by + 
for the upper Bide of the line segment and by — tor the lower aide. 



«?(*,) =0 

= f>H-- 
= 0 



c'<x,<l 



The surface potential associated with a vertical 
upward velocity of the air of magnitude w at the cle- 
ment located at x=Xi is 

^-2^^ (x-x l )*+( 3f +y0' dxi (5) 

where 

y=-Jl—x s and yi=-i/l—Xi* 

Equation (5) is fundamental to the description of the 
noncirculatory flow pattern since, by integration with 
respect to Xi, any admissible potential distribution may 
be obtained. The integrated result desired in this case 
as c'-»c is simply 



=vb&- log N(x,c) 



(6) 



where 



e—x 



Limit of (as C- c) 



c c 



-J c i 

Fiqubi 3— Representation of the sharp vertical displacement as a limit. 

The potential associated with bH is due to a vertical 
velocity distribution 

w=0 — l<ari<e 

= filb c<xKl 

and is (see <f> p> reference 1, p. 5) 



-X 



Xt— 1 



=-ftb^Ul =? cos-'c- (x-c) log N(x,c)] (7) 

Exactly similar considerations were made for the tab 
and lead to equations (4). 

It may be remarked that the analysis assumes no 
leak of fluid in the gap between the "aileron and the 
wing; that is, the gap is considered sealed. 

The following new sets of integral evaluations will be 
required. The expressions for the T and Y terms are 
listed in appendix B. The T terms are functions of c 
or of d only. When no explicit mention is made, e is 



X 
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M f =&£ (x-e)pdx 
This moment may be written 
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to be understood. The Y terms are functions of both 
c and d. 



4>nix ^(T s +2ljl=? cos-'e) 

fjtdx ^(T X -IT S ) 

j\(x-c)dx ^(Tt-lTo) 

jje( X -c)dx — ^czwro 

j%dx ^W-ITJ 

£4>,(x-c)dx= - h -f (Ts-lcjtt) 



Jf(x-c)dx "-f(T 7 -lTJ 

4>,dx ^(Y 1 +mT t ) 



X 

I 



£<f, y (x-c)dx ^ (Yi—mYs) 

j\(x-c)<fc ^ (Y t -mY t ) 

J]&& I^WlFr) 

f^-Z^-lYJ . 



X to{x-d)dx b ^(Y>-lY t ) 



x 

£<fc (x-d)dx= (Fa-iF,) 

The pressure difference on an element of the airfoil 
located at x is 



9 ( v i ^ 



(8) 



and the total force (positive downward) is therefore 

P=-2pbJ^4>dx (9) 
The moment on the airfoil (positive clockwise) about 



x—a is 



M a =6* J* ^ (x—a)pdx 

= — 2pb£ j>(x— c)dx+2pvbj ^<f>dx 

—2pb\e—a) j\ 4>dx 

The moment on the aileron (positive clockwise) about 
the hinge x=e is 



(10) 



M fl =— 2 P 6 l J 1 j>(x-c)dx+2pubj* <f>dx 

+2pb*l£ 4>dx+2ficUm\ (11) 

Similarly, for the moment on the tab about the hinge 
x=f, 

J a/ T =-2p6 , J^ 4>(x—d)dx+2pvb£ <f>dx 

+2pb*m£ 4>dx+2pcbm[<t>]l (12) 

Circulatory terms.— The potential at the surface of 
the airfoil associated with an element (counterclockwise) 
of vorticity in the wake at x 0 of magnitude AT= Ubdx a 
is (reference 1, p. 6) 



where x\<Cl and 2 0 >1. 
The potential for the entire wake is 



(13) 



(14) 



With the assumption that the wake remains where 
formed, the expression for the pressure at x (equation 
(8)) becomes 



V 2, 



6V cte 



Xo+X 



fdxo 



(15) 



i 2x71-^^-1 ° 
The Kutta condition for smooth flow at the trailing 

edge, which requires that d/dx($r+4>) must remain 

finite for x=l, leads to the result 



^J^U^j^^dx^ra+h+b^ay+^Tyo-lT^ 
+^(Tn-2lT l0 )bp+^(T 10 (d)-mT n (d))vv 

-r-^(r„(d)-2mr 10 (^)&7 (16) 

This quantity will be denoted by Q. 

When the various degrees of freedom of the airfoil 
are undergoing sinusoidal motions of the form for 
a long period of time, the wake is also sinusoidal and 
it is convenient to introduce the parameter k defined 
by the relation 



ut=ka=k^ 



or 



k= 



oh 



(17) 
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The force (negative lift) on the airfoil due to vortdcity 



is 



P=bppdx=-2xpvbO(jlc)Q (18) 

where C(k)=F(k)+iO{k) denotes the fundamental 
function introduced by Theodorsen (reference 1, p. 8). 
The moment about x—a is 

M a =b t j 1 _ i p(.x-a)cb=2TMf_(a+^C(k) -i]<2 (19) 
The hinge moment on the aileron about x=e is 



M„=6 2 ^ 1 2 >(J;-e)£iie=6 , V p{x-c)fe+b\c-e) C pdx 

Jc Jc Jc 

= ^pvb*(.T ls C(k)-T i )Q+pvbn(2T 2Q C(Jc) +2^/T^)Q 

(20) 

Similarly, the moment on the tab around its hinge, 
x=J, is 

+pcb*m (2r S0 (d) C(k) + 2 VT=^) Q (21 ) 

Finally, when both the noncirculatory and the 
circulatory terms are combined, there result the follow- 
ing expressions for the force and moments. 



TOTAL FORCE 

P=- p b\Th+VTa-*baa-vTj-T l bj3-vT < (d)j-T l (d)by) 

-pVmi^^T^Pvy+bTtid)!) -2*ptbCQ (22) 

where 

<2=m+A+&(i-a)d+£(:r 1( ,-mi)t>/j+^ 

TOTAL MOMENT ABOUT x=a 

M a =-pb s [-a*bh+v(Ji-a)vb<i+*b* (|+a J )a+T 1 ^i3+r ie t*a+2r v 6 t i 8+2 T «(^^Y+T , , s (^p67 

+2T li m i y^-pbn(T**(}+T n vbi}+T u b*(})- (23) 

TOTAL AILERON HINGE MOMENT 

\ Tir ir iTT ir/ 

-pb'lfTM+Tvvb a + T«W5 +-T 2t v l p+h r! bvH -T^H-Yn^y+^Y l2 bvy+-Y i b i y) 

-&pQTvi?p+±T*bv- l 3-lT i b*^-pb^ 

- f Anm(-Y 1 ^y+~Y 1 M--Y l b i y)-pvb 2 (.T lt -2lT il) )CQ (24) 

TOTAL TAB MOMENT 

M,= -p&f - T^&fc + :r 17 (d)cte+2T I3 (^ 

-h t ($b*j)-pb*m(T<(Qbfi+TK(.$vbi+Tu(<r)b 2 Ei+^ 

+ T T x (d)i?y+±T„Wbvi+lTMb i y)- pb^^T^^y+^Mbvy—T^b 1 ^ • 
- P b t l(-Y i ^p+-Y n bp$+-Y l b i p)- P bHm(-Y a ^+^ (25) 
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Discussion of the term T2&. — The concentrated sink- 
source representing the steep break (fig. 3) properly 
describes the main flow pattern, but the local flow 
pattern at the break is incorrect. The underlying 
theory excludes the possibility of representing the flow 
at a steep break. The limiting process therefore can- 
not be used in this simple theory, as far as the local flow 
is concerned. 

There is one term that depends on the local flow 
condition at the break. This term arises in the evalua- 

i-vi-4— r — < 

Figcee 4.— Representation of the effective mean camber line for a displaced aileron 
with binge position at x—r. as depending upon an additional parameter €. 

tion of [<f>p]l (equation (11)) and is present in the expres- 
sion for the force on the aileron. It occurs, then, in the 
expression for hinge moment Mg in the term T st 
associated with the coefficient Pfi. 

In order to picture the local flow and, at the same 
time, retain physical reality, it is necessary either to 
disregard a certain small neighborhood of the break or 
to spread the concentrated sink-source over a certain 
finite area. This end may be accomplished by regard- 
ing the mean camber line of the displaced aileron with 
rounded leading edge as depending on an additional 




.6 



Hinqe-offsef parameter, I = e-c 



parameter e' (fig. 4). Let e—c'=l' and c'— c=X. The 
velocity potential at the surface replacing equation (2) 
is then 

V 

where 4> e , e denotes ^[Vl— i 1 cos -1 c— (x—c) log N(x,ej] 

The analysis can be performed by use of this *equation 
and a similar one for <j>p instead of equations (2) and (3). 
The result will, however, differ essentially from that 
already presented in two respects: (1) The average value 
of I wfll be slightly less than e—c and will be nearly 

e—^-- (2) The term log N(c,c r ) will occur in T u 

to replace the infinite term log N(c,c). 

An effective v.alue of c' may be estimated in any 
given case or may be determined by experiment. The 
essential point is that the difference e' — c cannot 
become zero. It appears probable that c'—c is greater 
.than, say, 0.05 (e— c) and less than 0.40 (e—c). As an 
average value 0.25 (e—c) appears good. Within these 
limits the value of the term log iV(c,c') is fairly inde- 
pendent of the selection of c r . The effect of the choice 
of g' in the steady case on the hinge momentis illustrated 
in figure 5. Conversely, the experimental hinge- 
moment values may be used to estimate c'. 



.04 




O .05 .10 .15 .20 .25 .30 
Hinge-offsef parameter, 1= e-c 



(Steady ease.y-O.) 



(a) c-0. (b) c-0.5. 

Fkwbx 5.— Hinge-moment coefficient against hinge-offset parameter t-e— c for various values of the mean-camber-lina parameter X=e'- 
The ordinate la 

. -MtfTpVc'H-^Ta+TuiTri-H^Tn-iTaTn-TuTn) \-P^T»+2TnT„) 
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Equations (22) to (25) may be conveniently expressed 
in a coefficient form: 

P=- Trpv t V>(j ; A eh +aA e4 ,+fiA elt +yA cr y ) (22') 
M^-TpJb^Art+aAtv+pAa+yA^ (23') 



M fi =- Tp^b^A^+aA^ M«J+ 7^7) (24' 
M r =-* P <fb*(jA^+aA«,+[iA ti +yA i ^ (25' 



where 



w=kv(b, Ac^Ra+ilci, i4«=i?«+i/«. etc. 



and where 



Raa = 
Ra$= 
Rah = 
Ray = 

Rta = 
Rtf= 
Rth = 

R»y = 

Rea = 
Rtf= 

Rey 



-4„+^„+(|+«)[i;(2'„(<i)-2mI'„(<i))^'-J(r, D (i)-mr 1 ,(i))^] 
1 2G 

-B 7I+ ^_±(r 12 -2ir M )[i ? (r 11 (d)-2 m r I o(d))^ r -i(7' IO (d) - i»r„(d))^] 

~ /l \20.2F 
:-C n -±(T n -2lT w ) 2 ^+kT l0 -lT^ 



—O 20 
— G„- T 



Rda 

R<n= 

Riy- 

Iaa = 
I air 

I at 



— 2? rt -i(r u (d)-2inr J( ,(d))[(|-a)^-^] 
.= pD<*-^(.T„(a) -2m7' 20 (d))[^(7' Jl -2ir ll) )^ ! -^ {T 10 -lT n ) 2 ^ ] 



-z>„-^(r„(d)-2mr ai (d))^' 



-J\i+^*-£<rii(«9-2m7V,(d))^^ 

J{A«-(i+a)[i(T' 1 o-;r„)^ ! + 2 ^(r.i-2/r 1 o)2/]j 

^ r2 -(^+a)[^r 1 o(d)-m7' 21 (d))^+^ ;: (r tl (cO-2mr IO (d))2if']) 
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J -=iK+2^^- 2 ^[f r +G- tt > Jf ] 

/ w =^5«+^(r H -2ir 30 )[i(r m -;r J 0^+^(2 T n-2f7 T ia )2/|} 

*H[«rf¥t<j->F] 

/^=4[c'/ J2 +i(r 10 -Z7 i il )^+^(T 11 -2zr 1() )2i?'] 
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/. T =|[^+i(r I o(cO-mr 11 (^)^+^(r u ((0-2m7\ 0 (^)2 J F'] 
/ <fa =4jz? a2 +^;(2' I4 (d)-2m7' S o(d))[^+Q-a)2F]) 

/^=|{i>«+^(r I2 (^-2mr w (^)[i(r 10 -mi)^^(r I1 -2f!r 1 o)22f']| 



^(^(c*) -2mr„(d))2iP 



I*y=$D«+±(T a (<t) -2mr„(d))[i(r 1(( (d) -mr t i(d))^+^:(7 , n(d)-2mr t o(cO)2/J 



. ' 1 

^ 22*11 . 7?« 



Si=J(2 T 17 +Z2 I 2s ) 

706999 O - « - 16 



Bn=^{-T t +2lT t -VTA 

T 

S-ri=p(-F ( rKF 8 +mF t -Zml' 1 j 

Srt^CF^fF.+mFu+ZmF.s) 
CU=l 

c H =i(-i' t -2ivr : =^) 
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Dfi=B yl 

D 0 =p(Y a +mY ll +lY il +lmY u ) 

Z? Tl =i(- T s (d)+2mT 1 (d)-m i T li (d)) 

D*= ?(T K (d) +mT„(d)+m i T lt (d)) 
D^^^+mT^+m^d)) 



CONCLUDING REMARKS 

The material presented in the preceding pages repre- 
sents an extension of the work published in reference 1 , 
which has been expanded to include the tab functions 
and the effect of aerodynamic balance. Inasmuch as 
this addition fits in with the general arrangement of the 
earlier report, reference should be made to that report, 
and also to reference 4, for the application to the 
flutter problem. 



Lanqley Memohial Aeronautical Laboratory, 
National Advisory Committe fob Aeronautics, 
Langley Field, Va., December 16, 1941. 



APPENDIX A 
COMPARISON WITH REFERENCE 2 



In order to compare the results of this report with 
those of Kiissner and Schwarz (reference 2), the follow- 
ing relationships are noted: 



<■= — cos <p 
d=—cos 4> 
k=—io> 
C(k) = l+2T(-io>) 
1 

'2 



2i«-3»« 

-$17 



T s =- 



7i6=*i 



Tu 
T lt 
T u 

r„: 

T„ 

T, 

Tu 

T K 
T n 
T a 
T u 



1. 
t 

2 ' 



10 
= $18 

= -2* 5 



4*» 



T„= -(*,.+*») 
T*=i*r t * 



Tu=- 

T K =- 1 

log N(c,d) = ~Lb,t) 

r,M)=-.ar M (p,*) 

Y t {c,d) = Y lt (p&-T i {c)T n {& 

Y i (c/[) = -X l M^) 

Y t (e,d) = -X S M) 

r,(c,rf) = r,(cf,c) =X»(^) -Xvfaf) 

F T (c,oO +*u(*0*8i(*) 
F 17 (c,o0 = r,(o?,c) 
YM = Y l0 (d,c)=X,(<p,f) 

Y l9 (c,d) = Y ll (d,c) 

Y m (c,d) = Y lt (d,c) = -X<( v ,t) 
Y n (c,d) = Y a {d,c) = -Xnterf) 
Y a {c,d) = Y u (.d t c) = -X„(<p,+) 
Y»(c,d) = Y u (d,c) =X u ( Vl f) 
Y it (e,d) = Y lt (d,c) =X u ( l p,f) 

to be an error In the sign of the numerical value* for*u In table 3 



> There a; 
of reference 



APPENDIX B 

EXPRESSION FOB THE T AND T FUNCTIONS 



T FUNCTIONS 

T b =ct/ T=<* c os-'c- (l-c 2 ) 

=Vl-c* (c cos^c-Vl-c 2 ) 

2 , I =-i(2+c 1 )Vl :r ?+c cos" 1 c 

T 1 =c(l-c*)-(l+(?) : JT=c i cos" 1 c+c(cos- 1 c) 1 

7 T 3 =^-|(l-c J )(5c 1 +4)+ic(7+2c J )Vl ::i c i cos-'c 

-(l+c 1 ) (cos-'c)* 

z^eyi— c»— cos-v 

T s =-(l-G 2 )+2cVT^3 i cos" 1 c— (cos- l c)* 

r»= -|c(7+2c s ) JT=?-(±+c>) cos-'c 

r,= -|(1 +2c*) V W +c cos" 1 e i (l-c^-c^ 

r.=|[|(l-^ ,/l +a7;] 

2 T 10 =Vl zr c i +cos- 1 c 
2 , u=(2-c)VT ::: c i +(l-2c) cos^c 
J T u=(2+c)Vl = c 1 -a+2<;) cos" 1 c 

T u =~(T 7 +(c-a)T t ) 

r ls =r t +r u =(i+ c )VT=? 

r„= -2r,-r I +(a-i)2 T t = -c s ) ,/2 - r.-^r, 

T tt =T s -Tjr u 
1 

is — — 2 4-t II 

7 T »=— VT^+cos" 1 c 

r ss =2VT=?- A /^ 

r B =(— l-2c+2ffi)Vl = c i 

t m =t,+ (c-c)r«=-2r, 
r M =r 4 - (1-0)71-^ 



r iS =2Cl+c+log^(c > cQ) 
2jj=2-^l — cTib 

The term 7^ is discussed separately in the paper.. 
The variable c is to be understood when no explicit 
variable is indicated for the T terms. 

r FUNCTIONS 

F^^-Vl-eVl^-cos- 1 c cos" 1 d 

+d Vl-<? cos-'c+cVT^cos-'d- (<Z-e)*log JV(c,cO 
where 

1 -cd—y/T^-jL^d 1 



N(c,d) = 



~3=i 

=2-Jl—d 1 cos" 1 c— 2(d-c) log 7V(c,t2) 
r I (c,d)=i(c+2d5 VT^Vl^+rf cos^c cos-irf 
-"|(2+<P)yi^ cos" 1 c 

-la+SaZ-e 1 ) Videos" 1 d 

+i(c?-c) 8 logiv-(c,d) 
F 4 (c,d) = yi(rf,c) 

=|(d+2c)VT I= c s Vl I= ^+c cos" 1 c cos" 1 

— \{2+<*)->}T=<? cos-'c? 
-Ici+SctZ-d^Vl 11 ^ cos" 1 c 

-|(cf-c) 8 logiV( C> d) 

Y s (c,d) = —j\—<*ji—d*+ (2c-tOVT ZI ^ cos" 1 c 
+ (<*-c) , log2V(e,(*) 

— (g+^y cos -1 C COST 1 <£ 
+|[^(|-d l )+ C (2+d*)]vr= : ? cos" 1 c 
+|j(§-c i )+d(2+c 1 )]yi=? cos- 1 d 



+^ T ^logN(c,d) 
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2V 1 -c 2 cos" 1 d+2(d-e) log N(c,d) 

= Y i (d,c) 

Y s (c,d)=-- } /T^?^T^d i +(2d-cUl-e l cos" 1 d 

+ (d-c)>\og N(c,d) = Y,(d,c) 
Y*(p$ = Yi—T t (c) T w (d) 

Y w (c,d) = Y t - Y t -±T t (c)T u (d) 

Y u (e,d) = Y 7 - 2VT=?r 10 (d) 
Y lt (fi,d) = Fx- F 8 - Vl -^TnW 
F 1 ,(c ) d) = F 1 +2 T t (c)r 11 (d) 
F 14 (c,<Z)=F 5 -F 9 

F^fod) =2iJT=?T Si (d) +2 log N(e } d) 
Y lt (fi,d) = Y t - F T +2VT=?7 , I0 (^ 
F ir (c,d) = F,- T<(d) r 10 (c) = F 9 (d,c) 

F 18 (<!,d) = F,- F»- gl^I^e) =F I0 (d,c) 



Y lt (e,d)- 
F so (c^ = 
F 91 (<:,<*) = 
F M (e,d)= 
FaCe,^ 
Yu(e,d)-. 



= Y t -2iJT^T u = F u (rf,6) 

=F,- F.-VT^^i^ FhC^c) 

= F ? +7',(d)r 11 (c) = F Ii (d,c) 

= - F t + F„+ r 10 ( C ) r 4 (d) = F M (d,c) 

--2-JT=&T tl (c)+2 log N(fi,d)=Y ls (d,c) 

--Yi-Y t +2^¥T W = Y li (d,e) 



In the evaluation of the F terms the following per- 
tinent integrals occur: 

j(x-c)"logNdx=^^ \ogN 



where 



n+l J Vl 3 ? 



a;— c 



In order to evaluate the last integral put x=cos 6; for 
example, consider n—1: 



6) logNdz=-&j^- ] °g N ^ c > d ) 
£=?(JT=#-c cos- d) 
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